Constraining the nuclear equation of state at subsaturation densities 
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Only one third of the nucleons in 208 Pb occupy the saturation density area. Consequently nuclear 
observables related to average properties of nuclei, such as masses or radii, constrain the equation 
of state (EOS) not at saturation density but rather around the so-called crossing density, localised 
close to the mean value of the density of nuclei: p ~0.11 fm -3 . This provides an explanation for 
the empirical fact that several EOS quantities calculated with various functionals cross at a density 
significantly lower than the saturation one. The third derivative M of the energy at the crossing 
density is constrained by the giant monopole resonance (GMR) measurements in an isotopic chain 
rather than the incompressibility at saturation density. The GMR measurements provide M=1110 
± 70 MeV (6% uncertainty), whose extrapolation gives Koo=230 ± 40 MeV (17% uncertainty). 

PACS numbers: 21.10.Re, 21.65.-f, 21.60.Jz 



Constraining the nuclear equation of state (EOS) and 
reducing the uncertainties on nuclear matter properties 
in dense stellar objects such as neutron stars and super- 
novae is one of the major goals of nuclear physics. To 
do so observables such as masses, radii or energy cen- 
troid of the isoscalar Giant Monopole Resonance (GMR) 
are measured in finite nuclei. Relating the EOS to such 
observables is usually undertaken in several ways. The 
liquid drop expansion [l| is one of the most used methods. 
In its generalised version, one performs a development of 
a quantity (the mass for instance) considered at satura- 
tion density (volume term) , adding several terms such as 
the surface one. It should be noted that taking the vol- 
ume term at saturation density is based on the fact that 
the inner part of the nuclei density is close to the sat- 
uration. Another method is the so-called Local Density 
Approximation (LDA) [2|. In this approach the global 
properties in nuclei are obtained from considering nuclei 
as local pieces of nuclear matter. Finally, another method 
is based on the microscopic approach, relying on energy 
density functionals (EDF): using various EDF's, a corre- 
lation diagram is drawn between the predicted observable 
and a related property of the EOS. The measurement of 
the observable allows to validate an EDF and the corre- 
sponding property of the EOS For instance, the 
neutrons skin measurement is correlated with the slope 
of the symmetry energy. It should be noted that usually 
EDF are designed using masses, radii but also the nuclear 
incompressibility 0, Hi ■ 

Each of these methods comprise however limitations. 
The liquid drop expansion is known not to be accurate 
enough in the case of the incompressibility at the sat- 
uration density 0, [l(| , although the inclusion of higher 
other terms is relevant 11] . The liquid drop expansion 
of the mass has been successful, providing an accurate 
determination of the saturation energy. This quantity is 
now used in EDF determination, but this is an excep- 
tional case, related to the profusion of data available on 
masses. In the case of the LDA, its validity is generally 



questionable for finite size systems. Finally, in the case of 
the microscopic approach, the EDF employed has usually 
been adjusted to data on magic nuclei, whereas most ap- 
plications are aimed to be used for deformed open-shell 
nuclei. However there may be a general consensus that 
the microscopic approach should be used in fine, because 
of its accuracy and reliability. 

It seems for now difficult to straightforwardly deter- 
mine the nuclear incompressibility even with the micro- 
scopic method. The earliest microscopic analysis came 
to a value of K oo =210 MeV [|, but with the advent of 
microscopic relativistic approaches, a value of K oo =260 
MeV was obtained [l2j . It has been shown that this value 
is not determined accurately and that the density depen- 
dence of the EDF as well as pairing effects (and therefore 
the shell structure) have an impact on the determination 
of Kqo flljl - filil ]. Using Kqo in the design of EDF may 
not be a sound approach since it cannot be safely deter- 
mined by the microscopic method. Therefore the method 
of determination of the nuclear incompressibility has to 
be rethought and more generally, it is necessary to clarify 
the link between nuclear matter EOS determination and 
nuclear observables. 

First it should be noted that the liquid drop expan- 
sion is not a perturbative one since the surface properties 
of the nucleus are almost as important as the bulk one. 
Hence it may be a misleading idea to consider the nu- 
cleus as mainly composed of nucleons at saturation den- 
sity with a few at the surface. Let us consider the case of 
208 Pb which is usually considered as a benchmark for ex- 
tracting bulk properties. Fig. [T] shows the total density 
calculated in the Skyrme Hartree-Fock (HF) approach. 
The lower part shows the usual representation of the den- 
sity whereas the upper part displays an equivalent repre- 
sentation with an X axis scaled as r 3 . This allows to take 
into account the increase of nucleons per volume unit; 
the total number of nucleons corresponds to the integral 
of constant steps of the radial density represented on the 
upper part of Fig. [TJ It is now perceptible that about 
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FIG. 1. Matter density of 208 Pb calculated with the HF ap- 
proach using the SLy4 functional with different X axis scales. FIG. 2. (Color online) EOS incompressibility calculated with 
Top: representation taking into account the nucleons' distri- var i OU s functional, showing the crossing density around 0.7 
bution in the nucleus. Bottom: usual representation on ~ 11 fm -3 



one third of the nucleons of the 208 Pb nucleus lie in the 
saturation density area, whereas two thirds are localised 
in the surface at a density lower than the saturation one. 
Therefore even in heavy nuclei, the contribution of the 
surface is larger with respect to the volume one, raising 
the question of the legitimacy of constraining EOS quan- 
tities at saturation density by measurements of nuclear 
observables. 

Another illustration is given by calculating the mean 
density of the 208 Pb nucleus. Using a Skyrme-HF cal- 
culation, one obtains <p> = 0.12 fm -3 in 208 Pb with 
a variance \f <p 2 > — <p> 2 =0.04 fm -3 . Therefore the 
density value characterising a nucleus is not the satura- 
tion one (po=0.16 fm -3 ) but a significantly lower one, 
with a range spanning from po/2 to po- It should be 
noted that the dependence of the mean density on the 
nucleus is rather weak: 0.11 fm -3 for 120 Sn. In light 
nuclei, the mean density drops down as expected: 0.07 
fm- 3 for 40 Ca. 

Consequently the measurement of an averaged observ- 
able in a nucleus is more properly related to a correlated 
EOS quantity defined around the mean density than at 
the saturation density. This fact is illustrated on Fig. 
[21 where the density-dependent incompressibility, defined 
by[!,[ia 
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and obtained from various EDF's is plotted with sev- 
eral Skyrme, Gogny and relativistic interactions, all de- 
signed to describe observables in nuclei such as masses 
and radii. They intersect around the crossing density 
p c ~ 0.7po — 0.11 fm -3 . The existence of a crossing den- 
sity has been empirically noticed in previous works on the 



symmetry energy (Fig. 2 of Ref. [6]), pairing gap (Fie. 
2 of Ref. [3) or the neutron EOS (Fig. 2 of Ref. [|| 
and Fig. 1 of Ref. [H), and we provide here an explana- 
tion, related to the mean density: when designing EDF 
with nuclear observables, the corresponding EOS is con- 
strained not at the saturation density but rather around 
the mean density (the crossing density). In the case of 
the symmetry energy, the value at a density p ~0.11 
fm -3 is taken to be around 25 MeV, a value close to the 
symmetry energy coefficient of the liquid drop expansion 
[S II ' as an empirical prescription. This last value 
contains both a volume and surface terms, and thus rep- 
resents the symmetry energy extracted from nuclei ob- 
servables. For the incompressibility, the GMR is known 
to be related to the mean square radius of the nucleus by 
the energy weighted sumrule [2|. In the design of EDF's, 
the considered constraint on nuclear radii induces a con- 
straint on the compression mode, likely explaining the 
crossing around 0.7po observed on Fig. rjj This shows 
the universality of the crossing effect, arising from the 
constraints encoded in the EDF from nuclei observables. 
Due to this crossing area, a larger K oc =K(po) value for a 
given EDF can be compensated by lower values of K(p) 
at sub-crossing densities, so to predict a similar energy 
of the GMR in nuclei: the GMR centroid is related to 
the integral of K(p) over a large density range [17] . This 
allows to understand how EDF with different Kqo can 
predict a similar energy of the GMR, as noted in [l5[ . 

Various EDF's shall exhibit various density dependen- 
cies around the crossing point. At first order the deriva- 
tive of the pairing gap (or incompressibility or energy) 
at this point will differ between various EDF's. Comple- 
mentary measurements in nuclei are needed to charac- 
terise these derivatives. For instance in Ref. 0, Q the 
derivative of the neutron EOS around p c ~0.11 fm is 
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found to be constrained by neutron skin measurements in 
208 Pb. The associated interpretation was that the 0.11 
fm~ 3 density was considered as the neutron saturation 
density in nuclei. We provide an alternative and general 
view: 0.11 fm -3 corresponds to the crossing density, as 
seen on Fig. [2j 

We apply this method to the determination of the nu- 
clear incompressibility. When the Giant Monopole Reso- 
nance is measured and well reproduced by a given EDF, it 
shall therefore not be correlated with the incompressibil- 
ity of EOS at saturation density but rather with its first 
derivative M around the crossing density. It should be 
recalled that the crossing density exists because the EDF 
are determined by including nuclear observables such as 
masses and radii in their fit. On Fig. [2] the crossing 
point at p c ~ OJpo yields K(p c ) ~ 40 MeV for the in- 
compressibility (this is analogous to the fixed symmetry 
energy taken to be 25 MeV, as discussed above). To be 
consistent with a generalised liquid drop expansion [20| . 
the derivative M of the incompressibility at this point is 
defined by: 
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FIG. 3. Centroid of the GMR in 208 Pb and 120 Sn calculated 
with the CHF method vs. the value of M for various function- 
als. The experimental values for 208 Pb and 120 Sn are taken 
from Ref. [2J] and [2l[, respectively with respective error bars 
of ± 200 keV and ± 100 keV 



M = 3pK'(p)\ p = Pc (2) 

where K'(/j) is the derivative of the incompressibility 
density dependent term defined in Eq. (TTJ. 

Most of the GMR measurements have been performed 
on 208 Pb and data on other nuclei like the Sn isotopic 
chain is relevant [2l[ . In Fig. [3] the GMR prediction 
using the constrained Hartree-Fock (CHF) method as a 
function of M is shown for various Skyrme EDF's in 208 Pb 
and 120 Sn. The CHF method is a sum rule approach to 
calculate the centroid energy of the isoscalar GMR: 



-ElSGMR — . . (3) 
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The tyl\ moment is evaluated by the double commuta- 
tor using the Thouless theorem [2| and for the m_i mo- 
ment, the CHF approach is used 0, [H [Hj]: the CHF 
Hamiltonian is built by adding the constraint associated 
with the IS monopole operator. The CHF method has 
the advantage to very precisely predict the centroid of 
the GMR using the m_i sumrule. To be comprehensive, 
Fig. [3] also displays values for several relativistic and the 
Gogny D1S EDF's in the case of 208 Pb. 

The value of M deduced from both nuclei is compatible 
within few percents. In order to derive a sound value of 
M, we have further performed similar calculations in the 
ii2-i24g n i so t p eSi as we \\ as m 90 Zt and 144 Sm nuclei, 

leading to M~1100 MeV ± 70 MeV. The contributions 
to the uncertainty shall come from the small variance 
between the mean density in the respective nuclei com- 
pared to the crossing density, the isospin dependence of 
the incompressibility and the pairing effects. The linear 
correlation observed on Fig. [3] is striking, since the inter- 
actions employed can have a symmetry energy spanning 



from 30 to 37 MeV. This shows that the present results 
moderately depend on the isospin asymmetry of the sys- 
tem studied. The neutron vs. proton asymmetry pa- 
rameter <5=(N-Z)/A also remains rather constant in the 
considered nuclei for the GMR such as 208 Pb and sta- 
ble Tin isotopes: 6 ~ 0.2, validating the present isospin 
independent approach. 

It should be noted that the present method clarifies 
the issue of determining different incompressibility val- 
ues at saturation when using either the Pb or the Sn 
data [HI, [Hj]: a close M value is found with these two 
data sets. However the remaining discrepancy of the M 
values deduced from the Sn and the Pb measurements 
shows that the proper density dependence of the EDF's 
has not been revealed yet. Another striking feature is 
that the Gogny and relativistic EDF's are found on the 
same linear correlation than the Skyrme one, showing the 
universality of the ^gmr vs. M correlation, contrary to 
the ^gmr vs. one: it is well-known that relativistic 
EDF's can predict a similar Egmr in 208 Pb but with a 
significantly larger value of Koo jl2j ]. 

Let us now investigate why a clear correlation exists 
between the centroid of the GMR and in the specific 
case of the Skyrme functionals. [l5j]. Fig. 0] displays 
the relative contribution of the kinetic, central (to), fi- 
nite range (ti,t2) and density (t^) terms of the Skyrme 
functionals to the EOS and its derivative values at the 
crossing density. A striking regularity is observed among 
the functionals used. The dominant terms are the central 
and the density ones, but the central term vanishes from 
the second derivative of the EOS and beyond, allowing 
the density term (in p a ) to dominate alone. Therefore 
the second derivative terms and beyond are correlated 
to each other, implying that the correlation between M 
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FIG. 4. (Color online) Positive or negative relative contri- 
bution of the kinetic (circles), central (squares), finite range 
(diamonds) and density (triangles) terms to the equation of 
state (a), its first (b), second (c) and third (d) derivatives at 
the crossing density, calculated for 36 Skyrme functionals 



and Egma is propagated to the one between Koo and 
EGMi?: a linear correlation is preserved between M and 
due to the vanishing of the central term from the sec- 
ond derivative and beyond of the equation of state. In 
other words the density dependence of the incompress- 
ibility, driven by these derivatives, is correlated to M. 
Therefore is correlated to M, and this outcome is 
similar for all Skyrme functionals as seen on Fig. [4] the 
(Egm-RjKoc,) linear correlation may be due to an arte- 
fact from the specific density dependence of the Skyrme 
EDF's. 



But this linear correlation vanishes when including rel- 
ativistic EDF predictions, as mentioned above. M shall 
be the exclusive quantity which can be deduced from 
GMR measurements, whereas the determination re- 
lies on additional density dependence assumptions. Us- 
ing the M value at the crossing point, a linear expansion 
provides K oo =230 MeV. An uncertainty of ± 40 MeV can 
be inferred from the spreading of values on Fig. [2] 
obtained with the various functionals, which is of the or- 
der of ± 17%. It is therefore argued that measurements 
in nuclei constrain the EOS around the crossing density 
(namely the derivative of the EOS quantity at the cross- 
ing density) and deducing values at saturation density 
remains a mainly model-dependent extrapolation. In the 
case of nuclear incompressibility, it is proposed to change 
the usual Egmr vs. correlation plot (only valid in 
the specific density dependence of Skyrme EDF's) and 
to replace it by a more reliable and universal Egmr vs. 
M plot (Fig. [3]), where M is the derivative of the incom- 
pressibility at the crossing density. The measurement of 
the GMR in more neutron-rich nuclei [26|, H3] will cer- 
tainly open the possibility to study a part of the isospin 
dependence of the incompressibility. 
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